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The superfluid density near the superconducting transition is investigated in the presence of 
spatial inhomogeneity in the critical temperature. Disorder is accounted for by means of a random 
T c term in the conventional Ginzburg-Landau action for the superconducting order parameter. 
Focusing on the case where a low-density of randomly distributed planar defects are responsible for 
the variation of T c , we derive the lowest order correction to the superfluid density in powers of the 
defect concentration. The correction is calculated assuming a broad Gaussian distribution for the 
strengths of the defect potentials. Our results are in a qualitative agreement with the superfluid 
density measurements in the underdoped regime of high-quality YBCO crystals by Broun and co- 
workers. 

PACS numbers: 



I. INTRODUCTION 

The superconducting transition measured in real ma- 
terials is often smeared or broadened in temperature in 
a way that correlates with sample quality or disorder. 
A sharper transition is taken as a signature of a higher 
quality sample. On the other hand, a straightforward 
application of the Harris criterion^ implies that uncorre- 
cted disorder is irrelevant and does not affect the nature 
of the superconducting transition. The Harris criterion 
is modified for correlated disorder and implies that the 
transition can be broadened, depending on the nature or 
dimension of the correlation. 

In addition, to the situations considered by Harris, the 
existence of rare regions (analogous to Lifshitz tails 2 in 
the density of states of a disordered semiconductor) with 
higher than average critical temperatures may also sig- 
nificantly affect the properties close to T c leading to a 
smeared behavior of the order parameter as the transition 
is approached from the ordered phases. Such rare regions 
can occur with or without correlations of the disorder. 
This type of behavior is expected to be especially pro- 
nounced for cases in which the bare superconducting cor- 
relation length £o is short. High temperature supercon- 
ductors in the underdoped regime, having a short coher- 
ence length at low temperatures, present a good example 
of systems with increased sensitivity to various types of 
correlated disorder. Furthermore, many of these mate- 
rials exhibit twin boundaries, grain boundaries and/or 
disorder due to oxygen chains, all of which are examples 
of correlated disorder. Therefore, the high temperature 
superconductors present a good example of systems with 
increased sensitivity to various types of disorder. 

High-T c superconductors, because they are unconven- 
tional in the sense that the gap averages to zero over the 
Fermi surface, are very sensitive to disorder, although 
the sensitivity to uncorrelated random disorder is some- 
what mitigated by their short coherence lengths^ The 



role of disorder in superconductors has been an impor- 
tant subject of study for several decades^^. Imper- 
fections in the lattice structure are not only responsi- 
ble for the diffusive motion of charge carriers above T c , 
but also lead to non-uniformity of the attractive inter- 
action between them, ultimately giving rise to a spatial 
variation of the local critical temperature^' 7 . In addi- 
tion to structural irregularities, inhomogeneities in T c 
may be caused by the coexistence of superconductiv- 
ity and various density wave a 8 ' 9 ' 10 . Indeed, the simul- 
taneous treatment of several competing orders may be at 
the heart of a complete description of high-temperature 
superconductor a 1 1 i 12 and other strongly correlated sys- 
tems. 

Thus, obtaining a comprehensive understanding of all 
aspects of the complicated interplay between supercon- 
ductivity and various types of inhomogeneities, is defi- 
nitely an important goal. An experimentally relevant set 
of issues that one can specifically address in this context, 
would include the influence of inhomogeneities on the 
critical temperature, fluctuation corrections above the 
transition, as well as the behavior of the order parameter 
just below Tc- It is clear that if the characteristic length 
scale of the inhomogeneity is large compared to the T = 
superconducting correlation length, and the width of the 
distribution of critical temperatures is small compared 
to some average Tc, the problem can be studied within 
the framework of Ginzburg-Landau (GL) theory with a 
space-dependent critical temperature T c (r). The GL ac- 
tion is suitable for describing universal properties of a 
system in the vicinity of the critical point which are in- 
sensitive to the details of the microscopic Hamiltonian. 
However, the functional form of T c (r) is non- universal 
and is determined by the type of non-uniformity present 
in the system. 

In this paper, we study the behavior of the superfluid 
density near the transition in the presence of one spe- 
cific type of randomness in T c , caused by the presence 
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of widely scattered planar defects. The reason for this 
choice is two-fold. First, twinning planes, that may be 
present even in high quality crystals, can be regarded 
as planar defects locally increasing the critical temper- 
ature. If a twinning plane is located at x = 0, the 
Ginzburg-Landau (GL) action will contain a 5-function 
term —u6{x)\^f(T)\ 2 with u > O. 1 ^ Second, to study the 
effect of disorder in general, one needs to find stable so- 
lutions of the GL equation for arbitrary T c (r) and then 
average over realizations of disorder in expressions for 
the physical quantities of interest. This is difficult to 
accomplish because the GL equation is non-linear. The 
problem becomes more tractable, however, if the correc- 
tion to the GL action responsible for the change in T c 
has the form mentioned above, with u ^ 0. 

In this work, we consider a special case in which the 
spatial variation of the critical temperature is modeled 
by randomly located planar defects at points X%, each 
contributing a term 

u{v)\^>{v)\ 2 =~u6(x~X t )\y(r)\ 2 (1) 

to the GL action. Assuming that these imperfections are 
dilute, we obtain a correction to the uniform superfluid 
density that is proportional to the defect concentration 
rii . The calculations are performed at the simplest Gaus- 



sian level of GL theory. The disorder potential u is taken 
to have a Gaussian distribution with zero average value. 
The presence of defects that locally enhance the critical 
temperature, gives rise to a spatially decaying solution for 
the order parameter around them. The leads to a small 
but non-zero total transverse response at temperatures 
greater than T c , the critical temperature of the pure 
superconductor. The actual critical temperature corre- 
sponds to the region of strongest T c -enhancement and 
enters as a non- universal parameter in our treatment. 

Our work is motivated by recent experiments in the 
underdoped regime of high-quality YBCO crystals^ Ex- 
periments, done using the cavity perturbation method, 
show that, despite the high quality of the samples, ves- 
tiges of finite superfluid density persist even for temper- 
atures above ■ Our findings, based on the assump- 
tion of dilute planar inhomogeneities in the form of twin 
boundaries, qualitatively explain the results of measure- 
ments close to the critical temperature. 



II. METHOD OF EXPANSION IN POWERS OF 
THE DEFECT CONCENTRATION 



The starting point of our calculations close to criticality is the Ginzburg-Landau (GL) free energy in terms of the 
local superconducting order parameter ^(r)*^ 

T = T n + J dr|^-|V*(r)| 2 +a(T)|vI/(r)| 2 + C/(r)|*(r)| 2 + ^(r)| 4 |. (2) 



Here !F n is the free energy of the normal system and 



a(T) = a 



T-T, 



(or 



(0) 



(3) 



is the distance from the critical temperature, Tj ', of a 
homogeneous superconductor. We assume that the de- 
viations of T c (r) from Tc°\ described by U(r), occur in 
regions of a size greater than or of order the T = corre- 
lation length, £o = £(T = 0), but small compared to the 
correlation length near Tc°\ This assumption justifies 
the use of the GL formalism for both conventional and 
short coherence length superconductors, provided that 
they are not too far from the critical temperature. In this 
case, if the regions where the critical temperature differs 
sufficiently from are located around points Ri, we 
can quite generally model the randomness in Eq. (|2|) by 



U(r) = $^u i (r-R i ). 



(4) 



In the subsequent treatment, we will refer to these re- 
gions as defects. The functions ttj(r — Rj), that we will 



henceforth call the potentials, are presumed to be quickly 
decaying with |r — FLj|. From the above considerations, 
the characteristic lengths of this decay must exceed £o 
but be small compared to the correlation length £(T) 
close enough to the critical temperature. In addition, 
it must be much smaller than the average separation be- 



tween the positions of the defects n 



-l/d 



In the Gaussian 



approximation of GL theory f(T) = 1/ y / 2ma(T). 

The equilibrium distribution of the superconducting 
order parameter, \&o( r )j follows from the solution of the 
saddle point GL equation that is derived by varying 
Eq. |0 with respect to **(r): 



-Yl + a(T) + U(r) 
Zm 



tfo(r) + b|*o(r)r*o(r) =0.(5) 



For a given disorder potential, the actual transition tem- 
perature, T c , is determined from the value of a(T c ) for 
which a non-zero solution of \l/o(r) first appears. This 
happens when the eigenvalue spectrum of the operator 



Lo[U] 



2m 



+ a{T) + U(r) 



(6) 



3 



crosses zero. The ensuing distribution Vt^r) can be cho- 
sen real and positive everywhere in space and must be 
stable. The stability conditions can be determined if 
one expands the generally complex order parameter ^(r) 
around *o(r) in Eq. ©, 

*(r) = * (r) + *y(r) +i*x(r). 

One can then easily infer that a non-negative eigenspec- 
trum for the operators 

Lx[U] = -Z-+<x(T) + U(r)+b* 2 (r), (7) 



L l[ [U] = -—+a(T) + U(t) + 3b^l(v), (8) 

is a necessary condition for stability. The problem of 
determining the transition point and finding the stable 
solution below the transition analytically for the general 
form, U (r), is a daunting task. However, one can simplify 
the problem if the concentration of defects n, is small. If 
one assumes that every defect potential, «j(r — Rj), is 
characterized by the same set of parameters {it}, one 



can employ the method of expansion in powers of con- 
centration riit^ 

Consider the function 

Fjv({ui}Ri, • ■ • , {wat}Rat; r) 

that describes the spatial dependence of some quantity of 
interest and is calculable based on the GL action in the 
presence of N defects. We presume also that every defect 
located at point R^ has a potential characterized by the 
specific parameter set {ui}. The values of the parameters 
{u} are distributed according to the distribution V({u}), 
satisfying the normalization condition 

J V({u})V{u} = 1. (9) 

In Eq. formally means the integration over all 

variables in the set {u}. We are interested in the value 
of the function F/v({mi}Ri, . . . , {wjv}Ra?; r) that is an 
average over positions Ri as well as parameters {ui}. If 
one denotes the positional average by angular brackets, 
the full average can be written in the form: 



F N ({u 1 }R 1 ...{u N }R N ;r) = J V({m}) . . . V{{u N })V{ux] . . . V {u N }(F N ({ui}Ki, . . . , {u N }R N ; r)>. (10) 

Regarding the concentration of defects as small, one can formally write the average Eq.fTOf in the form of a series in 
powers of rij. The details of the corresponding derivation are presented in Appendix A. As follows from Eqs. (|A2[) - 
(|A4[) , up to the first order in rn, 

(F N ({u 1 }R u ...,{u N }R N ;r)) = F (r) + n i J [Fi({ui};Ri; r) - F (r)] dR u (11) 

so that 

F N {{u 1 }R 1 ...{u N }R N ;r)=F (r) + n l J V({u})V{u} J [Fi({u}\ R; r) - F (r)} dR (12) 

Eq. (fl2|) contains the lowest order correction to the function Fo(r), the quantity of interest in the absence of any 
defects. The calculation of this correction requires the knowledge of function Fi({u}; R; r) - the quantity of interest 
in the presence of just one defect located at point R. We should mention that the approach described in Appendix A, 
provides a way to reduce the level of complexity of the initial problem, since it reduces to calculations in the presence 
of just a finite number of defects. This task is simpler, although in practice one has to limit the treatment to the 
level of one or at most two lowest orders in m- An important assumption made in the development of this approach 
is that all integrations in Eq. (| A4|l . in every term of expansion in m, do not lead to divergences. Convergence must 
be maintained for all parameters in the set {u} and all values of ct(T) especially the point a(T) = 0. This property, 
ensuring that this method of expansion is controlled, is far from being a forgone conclusion, and must be carefully 
addressed once the specific form of the defect potential is chosen. As will be shown below, for the essentially one- 
dimensional potentials such as those given by Eq. ([1]), the convergence is maintained for all a(T) in the first order 
of expansion in rij. Although, we will be concerned below only with this lowest order, we believe that the procedure 
is well-behaved at all orders, as long as defects are parallel to each other and the one-dimensional character of the 
problem is maintained. 
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III. SUPERFLUID DENSITY 



In this Section, we calculate the superfluid density in the presence of randomly located planar defects based on the 
expansion in powers of defect concentration n.j. We will limit ourselves to calculations up to first order in n.j. It is 
useful, however, to first discuss the qualitative behavior of the superconducting order parameter in the presence of 
defects, without specifying the dimensionality of the problem or the form of the potential u(r — Rj). In the absence 
of randomness, the superfluid density isi^ 



pf(T) = [^(T)f = | w - |/6 



0, 



a(T) > 
a{T) < 0; 



(13) 



To calculate the first order correction in rii to this result, one needs to solve the saddle point equation for the order 
parameter ^/^(R;; r) in the presence of one defect located at point R ; : 



■— +a(T)+u(r-R l ) 
2m 



^ 1) (R i ;r)+6 *^(R<;r) *^(R i; r) = 0. 



v(D 



(14) 



It is clear that (Rj; r) = (r — Rj), and without loss of generality we can consider the defect to be located at 
R s : = 0. The non-zero real and positive solution of this equation occurs at the point where the eigenvalue spectrum 
of the operator LqM, containing a single-defect potential, reaches zero. In analogy with Eq. ([6]), 



L [u] = -—+a(T)+u(r). 
Am 



(15) 



Provided that the solution obtained from Eq. (|14|) is stable, we can write down the general formula for the superfluid 
density averaged over the randomness in u(r) 



Ts{T) 



V({u})V{u} l[^\r)] 2 dr, 



a(T) > 0; 



\a(T)\/b + ni JV{{u})V{u) j{[^\v)f-\a{T)\/b}dv, a(T) < 0. 



(16) 



r 



Expression Eq. (|16p follows straightforwardly from 
Eq. (fT2|) . applied to the square of the order param- 
eter. But the same result can be obtained using 
time-dependent Ginzburg-Landau equation as a starting 
points One needs to calculate the transverse response, 
and subsequently perform the average using the same 
method of expansion in powers of rii in the long wave- 
length limits The meaning of notations J c and V({u}) 
in Eq. ([16)) will be discussed below. We only mention now 
that we must carefully integrate, not over all possible val- 
ues of parameters from the set {«}, but only over those 
realizations that, first, give stable solutions for ^^(r) 
and, second, result in a defect-affected critical tempera- 
ture not exceeding some fixed value T c . 

If the system contains defects that give rise to a sta- 
ble positive non-zero solution of Eq. ([M]) at T greater 

than Tc°\ superconductivity must be presumed shifted 
to higher temperatures. In this case, a finite transverse 
response will be observed above the critical temperature 
of the homogeneous sample. To first order in rii, the ac- 
tual transition point will be determined by the defect that 
leads to the strongest enhancement of T c in the sample. 
This means that in experiments, the temperature T c , at 
which the superconducting response is first seen will be 



disorder-dependent and non-universal. Very close to T c , 
the superfluid density will be tiny, since the order param- 
eter will be determined by contributions coming from a 
very small number of defects. But with decreasing tem- 
perature, the fraction of defects giving rise to non-zero 
solutions of Eq. (Ti"4"|) will increase, leading to an increase 

in the superconducting response. At T > T c (0) , the so- 
lutions for the order parameter, ^^(r — Rj), will be 
localized around the center of the defect at point Ri. In- 
deed, as long as a = a(T) > 0, in the absence of any 
defects, the only stable solution for the order parameter 
is zero. Hence, it follows from Eq. (fT4")) that for functions 
u(r) that vanish quickly enough with distance, 

*W( r ) = ^£) e -v^I^ (17) 
Vb 

when the condition *Jmur 3> 1 is satisfied. The func- 
tion <7>(r) which has a weaker than exponential depen- 
dence on r, is determined by the effective dimension- 
ality of the problem. For instance, in two dimensions 
g>(r) = 1/y/r and does not contain any variables re- 
lated to the potential^ For ^fmar < 1, however, the 
functional form of the decay of the order parameter is 
no longer exponential and depends strongly on the non- 
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universal characteristics of the function u(r). Eq. (fT7|) 
also describes the behavior of 'J/q (r) at all distances 
when a — exactly. The corresponding asymptotic 
forms are in fact written out in Table 1 of Ref. LL8| and 
we will not discuss them further here. 

When T < Tc and a < 0, the stable solution for 
the clean system is given by ^^(T) = yj\a\jb, and we 
expect from Eq. (TTi)) . that finite values of u(r) will add 
some perturbation to this solution that falls off at infinity. 
We look for solutions of the form 



*^(r) 



1 

7h 



|a|+V(r) 



(18) 



where the real auxiliary function ijj(r) satisfies the equa- 
tion 



2m 



+ 2\a\ +u(r) 



ij)(r) + y/\a\u(r) 



-3V\a\mr)} 2 + [^(r)] 3 = 0. 



(19) 



By analogy to Eq. |[PT|) and provided that yjm\a\r ^> 1, 
we can write 



V (1) (r) =p<(r)e- 



(20) 



with the function g < (r) having an asymptotic form sim- 
ilar to that of <7>(r). Again, closer to the defect when 
y/m\a\r ~ 1, the crossover to a different functional form, 
with stronger dependence on characteristics of the po- 
tential, will take place. Note also that the requirement 
for V?q (r) to be positive, does not prevent the function 
g<(r) from having both signs. This means that depend- 
ing on the form and sign of u(r), the order parameter 
may be either enhanced or suppressed in the vicinity of a 
defect for a < 0. Thus we conclude that the spatial vari- 
ation of the order parameter changes qualitatively when 
a passes through zero, meaning that the defects play a 
different role in the system above and below . Above 
Tc°\ only a portion of all defects will perturb the zero 
value of and, though widely scattered, they nev- 
ertheless ensure a small but finite superconducting re- 
sponse. At the same time, below Tc°\ every defect will 
affect the solution ^ Q ^ = y/\a\/b, but this just leads 
to a correction to the superfluid density that becomes 
more and more innocuous with decreasing temperature. 
In some sense, we can say that, because of these qualita- 
tive differences, the point T = acts as a special kind 
of critical point. Indeed, a simple inspection of Eq. (fl6|) 
reveals that, despite the continuity of the superfluid den- 
sity at a = 0, its derivative with respect to temperature 
exhibits a jump. 



It is then appropriate to ask what kind of a defect po- 
tential u(r — Rj) leads to an increase in the critical tem- 
perature. To answer this question, consider the eigen- 
value problem for the operator Lq[u] — a, where Lq[u] is 
given by Eq. ([±5]) : 



V 2 , 
-— +u(r- 
2 m 



R, 



Xe(r-R 4 ) = e Xe {r-R t ). (21) 



This equation is nothing other than the Schrodinger 
equation determining the stationary states of a particle 
moving in the presence of potential u(r — Rj). For a po- 
tential which falls off rapidly enough at infinity, the spec- 
trum of positive eigenvalues e is continuous. It is describ- 
able by a number of quantum variables, with k — v2me 
being one of them. The spectrum of negative eigenval- 
ues E n , if any exist, is discrete. All eigenvalues e are 
explicit functions of all parameters in the set {u}. It is 
easy to see then that the transition point, in the pres- 
ence of one defect, is determined by the smallest eigen- 
value, Eq, of the operator Lq[u] — a. From the condition 
a — \Eq\ = and Eq. © it follows that, because of the 

defect, T c — Tc (1 + \E \/a). If there are no discrete 
levels, Eq — and no increase of the critical temperature 
occurs. We note that the points of instability coincide 
with the poles of the Green's function of the operator 
Lq [u] , obeying the following equation 



Lo[M]0W(R i ;r,r')=<5(r-r'). 



(22) 



The superscript means that the Green's function are cal- 
culated in the presence of only one defect located at Ri. 
Above T c , 



(Rj; r, r') = V X«(r - RQxjfr' - RQ 



(23) 



where S{e} formally denotes the summation and inte- 
gration over the discrete and continuous branches of the 
spectrum. Since any defect breaks translational invari- 
ance, (R,; r, r') depends separately on r and r' rather 
than on r - r'. The Green's function in the momentum 
representation (Rj; p, p') is the Fourier transform of 
Eq. (|2"3"|) . It contains two momenta, p and p', and has 
the same pole structure but may be more straightforward 
to calculate depending on the form of potential. 

Following the same strategy that led to Eq. (TTS)) . it 
is possible to write down the Green's function averaged 
over the positions of defects and parameters of their po- 
tentials. In analogy with Eq. (Tl2^) . up to first order in 
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g( N H{u 1 }R 1 ...{u N }R N ;r,r>)=gW(r-r>) + n l Jv({u})V{u} J |^/ (1) ({u}; R; r, r') - £ (0) (r - r')J dR (24) 



r 



An approximation that contains only the first power in 
rii is known as a single-site approximatio n 19 ' 20 ! 21 . If all 
defects have one and the same potential, the system is 
considered to have binary disorder, in which case the in- 
tegration over T>{u} would be absent in Eqs. (fT2"|) . (l2"4")) . 
For this widely studied simplified type of randomness, 
the Green's function Eq. (f2"4")l is obtainable as a result of 
resummation of a certain class of diagram o 19 ' 22 . If the 
distribution of parameters characterizing the potentials 
u(r — Rj) is broad enough, averaging over them with the 
weight V({u}) plays an important role introducing an 
additional complicating ingredient to the problem. The 
critical value of ot{T) for which the non-zero solution of 
Eq. (|14[) first appears in this case, will properly coin- 
cide with the singularity in the disorder averaged Green's 
function as can be seen from Eqs. (|2"5|) , (j2~4"|) . It will be 
determined by the defect that induces the maximum local 
T c in the sample. 

Next we address the question of the stability of possi- 
ble solutions of Eq. ([T4]) . In the presence of one defect 
(presumed located at R, = 0), the solution is stable if the 
eigenspectrum of operators and £j_ and ey, 

written in analogy with Eqs. ([7]) and {§]), is non-negative. 
That is, one needs to analyze two equations: 



~ + a+ u(r) + b[*P(T)]*\ f ± (r) = s ± f ± (r), (25) 
2m 



2m 



+ a + u(r) + 36['4 1) (r)] 2 (r) = en/u (r). (26) 



To do this, we employ the following general mathematical 
result for the spectra of second-order differential opera- 
tors. The eigenvalues can be ordered in a sequence of 
increasing values, and the eigenfunction corresponding 
to the lowest eigenvalue (the ground state eigenfunction) 
has no nodes as a function of Eigenfunctions corre- 
sponding to higher energies must change sign somewhere 
in space and are orthogonal to the ground state eigen- 
function. Comparing Eqs. (fl"4"|) and P5|) . we see that 
the eigenfunction corresponding to e± = is given by 
^/^(r). Hence if it is everywhere positive, we can claim 
that it is the ground state of operator L±[u}. This re- 
sult immediately implies the conclusion that the lowest 
eigenvalue of Eq. (j2"6")l , £m , cannot be negative. Indeed, 
the ground state eigenfunction /i|o(r) must be bounded 
and can not change sign anywhere. Hence if we consider 
the ground states of Eqs. ([23]) and (j2l)|) , multiply them 

respectively by /||o(r) and \Dq (r), integrate over dr and 
then subtract the first from the second, we find that 



2b 



<4 1} (r) /|| (r)dr= e 



II o 



This immediately implies that e» > is the only pos- 
sibility, and hence that all other em > as well. This 
result is just a simple manifestation of the fact that, for 
predominantly positive potentials, the set of eigenvalues 
shifts up. We can not say, however, whether the ground 
state belongs to the discrete spectrum or lies at the bot- 
tom edge of the continuous one. From Eq. (f2"6"]l it fol- 
lows that for a > 0, the continuous spectrum starts at 
En = a, while for a < 0, £|| = 2\a\ is its lowest possi- 
ble eigenvalue. If < eyo < a and a > (or 2\a\ for 
a < ), then the ground state eigenfunction /||o(r) be- 
longs to the discrete branch of the spectrum and falls off 
exponentially at infinity. But if emq is the lowest possi- 
ble eigenvalue of the continuous spectrum, /yo( r ) tends 
to some non-zero constant as |r| — > oo. Similar analy- 
sis of asymptotics following from Eq. (|25[) leads to the 
result that for a < 0, the spectrum of L±[u] is purely 
continuous and starts right from the ground state zeroth 
eigenvalue. If a > 0, however, one can not exclude the 
presence of some additional energy levels belonging to 
the discrete spectrum in the segment < s± < a. To 
conclude, if one finds a solution of Eq. (fl4"|) &q (r) pos- 
itive for all r, it is guaranteed to be stable. We are not 
aware of any general analytic methods that allow us to 
solve Eq. (fT4")) because of the cubic non-linearity. But 
the possibility that solutions may be found, depending 
on the relation between a and the parameters of it(r), 
seems quite realistic for smooth potentials of a given sign 
which decay monotonically at infinity. 

We are now in the position to discuss how to perform 
the average over disorder realizations J c in Eq. (|16[) , and 
what limitations one should impose on the distribution of 
randomness in order to obtain physically sensible results 
for the superfluid density. As discussed earlier, our ap- 
proach implies that the actual critical temperature is de- 
termined by the defect which gives the greatest increase 
of ri°^ . This approximation is a consequence of consid- 
ering the problem in the lowest order in nt and suggests 
that, to compare the theoretically calculated disorder- 
smeared behavior with experimental data, the actual T c 
should be introduced by hand. This can be achieved if 
the integral J T>{u} is performed, not over all possible 
values from the set {u}, but rather over those of them 
that do not allow for the defect-shifted critical temper- 
atures greater than the stipulated T c . Subscript C in 
Eq. (fT6|) is used to indicate exactly that. Since in this 
case J c V({u})T>{u} < 1, it is appropriate to introduce 
the normalized distribution, 



V({u}) = ^V({u}), A = jv{{u})V{u} (27) 
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and employ it in the average over disorder potentials. 

In addition, to eliminate the possibility of unphysical 
behavior of the superfluid density as a function of tem- 
perature, the disorder distribution must be regarded as 
symmetric about its average. Namely, we must require 
the fulfillment of the constraint: 



T({u})T>{u} / u(r)rfr = 0. 



(28) 



To clarify its importance, we first note that the positive 
sign of pi does not follow automatically from Eq. (jT4j) 
when a < 0. Indeed, if we divide Eq. (| 14[) by 6\&q (r — 
Hi) and integrate over space and T>{u}, then using Eq. 
(Unj) for a < 0, we find that 



Ps b 



^ jv{{u})V{u) J 



V 2 ^ x) (r) 
_2m^\v) 



u{v) 



dr (29) 



We can then integrate by parts the term containing 
^'(r), with the help of Eqs. (JTHJ), CP- Hence, if 
Eq. ([28]) is satisfied, it follows that 



|£V| 

b 



2mb 



V({u})V{u} 



V^ 1} (r) 
*«(r) 



dr (30) 



that pi > everywhere below . This could not be 
the case, had the left hand side of Eq. ([25)1 been negative. 
Although nothing wrong occurs if it is positive, it is con- 
venient to ensure that Eq. (|2"8|) is satisfied by adjusting 

Tc°^ which so far has been assumed to be the transition 
temperature of a disorder-free sample. 

We now apply this general formalism to the case in 
which the defects are described by the potential: 



u(r) = — u5(x — Xi) 



(31) 



The integrand in Eq. ([30]) is always positive, ensuring 



As mentioned in the Introduction, our model consists of 
a stack of parallel planes that locally change the critical 
temperature. The planes are infinite in y and z direc- 
tions, resulting in an essentially a one-dimensional prob- 
lem. Eq. (|3~lj) also implies that the planes are formally 
of zero thickness. Physically, this corresponds to a situa- 
tion in which the actual width of the planar defects is of 
order £o- The only new parameter in the problem having 
dimensions of length is l/(m|u|) which, together with the 
correlation length £ = l/v2ma, determines the charac- 
ter of the solution for the order parameter. Comparing 
those two length scales one can see that any shift to a 
should be proportional to mu 2 . The order parameter 
^0 ( r ) depends only on x, and assuming, again without 
loss of generality, that the defect is at the origin, we write 
the equation 



2m dx 2 



+ a — u5{x) 



= o. 



(32) 



The corresponding solution (x) must be continuous, but its first derivative has a jump at x = 0, meaning that 



dx 



z=+0 



dx 



-2mKf^(0) 



(33) 



x=— 



As has been discussed above, two qualitatively different solutions are possible depending on the sign of a, and we 
must consider separately two cases. 



A. Case a > 



where the constant Ai , determined from the condition of 
Eq. ([33]). is given by 



In this case, the solution and its first derivative must 
decrease exponentially at infinity. Elementary integra- 
tion then leads to the result that 22 24 



1 2a 

Ai = arctanh— \ — , (35) 
u \ m 



Note that the solution Eq. ([M]) makes sense only if u > 0. 
At large |x|, 



I2a 



\/fesinh [\/2ma|x| + Ai] 



(34) 



*W(x) 



f8a ( Vr 

T [7; 



/2a 



1/2 



/2ma|a:| 



2a 



, (36) 
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in agreement with Eq. (|17p. while at 




FIG. 1: (Color online). Order parameter v&q 1 ' plotted as a 
function of x using Eq. (|34[1 for m = 2.0, b — 1.0 and u — 0.5 
and two values of a shown in the figure. 



2a 



(37) 



This solution is possible only if < 2a < mu 2 . If 
2a > mu 2 , the only stable solution is (a;) — 0. We 
conclude that the amount by which u increases the crit- 
ical temperature is connected to a by 



(38) 



It is instructive to check this result by calculating ex- 
plicitly the Green's function form Eq. (|2"2")1 and finding its 
poles. It is convenient to work in the momentum repre- 
sentation. Since the problem is translationally invariant 
in the y and z directions, we seek solutions of the form 



j POO />00 



dq y dq z 



OO J — OO 



OO />0O 



ib*+pV) gU(r) yz ;p,p')dpdp', 



(39) 



— oo J —oc 



where r\ yz = r)(q y ,q z ) = q 2 /(2m) + q 2 /(2m). Substituting into Eq. (|22[) and performing simple manipulations in the 
term involving the (5-function, we obtain two equations to be solved self-consistently 



The result for the full Green's function is 

,(!)/„ ,s $(p + p') u 



C( P ') 



n ,,dp" 



G ( '(Vyz]P,p) = 



1 



2vr 



a + i] yz + p 2 /(2m) 2tt[1 - uK] [a + ri yz + p 2 /{2m)] [a + rj yz + (p / ) 2 /(2m)} ' 

I 



(40) 



(41) 



in terms of the integral 



K = K(j], 



dp" 1 

~2^ ' a + Vyz + (p") 2 /(2m) 
m 



v /2m(c 



Vyz) 



(42) 



B. Case a < 

In this case, as \x\ — > oo, (x) — » y/\a\/b. It is easy 
to complete the integration to obtai n 22 ! 24 



•\/|a|/6coth 



-\/|a|/6tanh 



A/m|a||x| + A 2 



\/mH M + A 2 



, u > 0, 

( 

, u < 0, 



(43) 



The first term in Eq. (|41[) is the Green's function corre- 
sponding to the absence of any potential, and has the 
simple pole at a = 0. The second term is the non- 
translationally invariant contribution due to the presence 
of the defect. For positive u only, the factor containing 
K has an additional pole at a = mu 2 /2, implying an 
increase of the critical temperature in agreement with 
Eq. 



. 1 . , 1 4\a\ 

X2 = -arcsmh: — -\ . 

2 \u V m 



(44) 



For a < 0, potentials with both signs of u lead to 
p hysica lly sensible positive solutions. At large distances, 



mu ■ e 



-2-v/ m|a| |a;| 



/ \a~\ + y/\a\ + mu 2 /A 



> (45) 
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FIG. 2: (Color online). Order parameter plotted as 

a function of x using Eq. (|43[1 for m — 2.0, b = 1.0 and 
it = 0.5 (upper figure), u — —0.5 (lower figure). The values of 
a corresponding to each curve are displayed in the figures. 



which has the asymptotic form discussed in the previous 
Section. Full expressions for * 1} (0) can be derived in 
a straight forward manner, but here we only present the 
less cumbersome ones in the limit of large and small (with 
respect to mu ) \a\, 



m a 



*W(0) 



u < 0, 



mu 2 > \a\. (47) 



u > 0, 



Depending on the sign of u, the order parameter is ei- 
ther greater or smaller than \E , g°' ) (x) = y/\a\/b. Thus, if 
a tends to zero from below, the solution for the order 
parameter vanishes if u < 0, but transforms into 



0. 



(48) 



b mu\x\ + 1 

for positive u. The dependence on x becomes a power 
law, meaning that the influence of the defect is long range 
at the special point a = 0. 

We can now substitute the solutions given by Eqs. (|34|) . 
(j43| into Eq. (|16|) and calculate the average superfluid 
density pl(T) . u is taken to obey the symmetric Gaussian 
distribution, 



P[u] 



1 



2ttW 



exp 



u 2 \ 
2W 2 J ' 



(49) 



with mean W. This form of distribution implies an 
exponentially rare probability of occurrence for defects 
with potentials with strength much greater than average. 
However, since u can, in principle, take any value in the 
distribution Eq. (|49|) , the defect-induced enhancement of 
the critical temperature is formally unbounded. We must 
then impose an upper limit u m on possible values of u, 
which will define the actual critical temperature 



T c = T r (0) ( 1 + 



2a 



(50) 



*^(0) 




t a\ I u I m 

T l 1+ 2V Id 



The renormalized distribution to be used in Eq. (|16p is 

P[u}du, (51) 
mu 2 \a\; (46) and it is simple to integrate over x and obtain 



pJCO = 



|0| 

b 



2rij 
b 

2n,j 
~b~ 



P[u] 
A 




P[u] (cothA 2 -l) 2 



A 



coth A2 



du, 



a > 0, 



a < 0. 



(52) 



In deriving the result for a < 0, we explicitly used the symmetry, P[— u] = P[u]. Noting that, from Eq. ([44]) . 

coth A 2 = ( Vjjl^+^ j) V2 , (53 ) 
V y4|a| + mu 2 — i/m\u\ J 
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we see that for a = the two expressions in Eq. (|52|) 
become identical. The presence of the Ti°^-enhancing 
defects makes the superfiuid density finite at that point. 
Estimating its order of magnitude, we can write 



mW 



(54) 



The value in the righthand side of Eq. (|54|) contains the 
first power of rij and is presumed small enough so that 
the whole approach based on the Ginzburg-Landau ex- 
pansion remains valid. 

For convenience, we introduce the rescaled parameters 



m 

2^ 



— u — > u, 
2 



W, \ —m 

m 



and plot 6pj(T) as a function of a(T) = a{T-T^ ] )/T^\ 
assuming a = 1 . The results for several values of disorder 
distribution width W ( W = 0.0, W = 0.1, and W = 0.2 
) and u m = on all plots, are presented in Fig. [3] 

From the plots it follows that, if the values of u are 
broadly distributed and u m ^> W belongs to the Lifshitz 
tail, the superfiuid density for larger a is exponentially 
small. Under the same circumstances, the behavior near 
the point T = is not sensitive to the exact value 
of u m . We note also that for u m 3> W, A f» 1 and 
is not of much importance. Below Tc , however, the 
behavior asymptotically approaches that of the disorder- 
free system. These results are in qualitative agreement 
with the solution for the order parameter presented in 
Ref. [13, having a small but finite value going to zero at 
some weakly ^-dependent value T c > . 

In Fig. SI we present data for experimental measure- 
ments of the superfiuid density taken from Ref. 0- The 
upper figure shows the results for the superfiuid density 
of an YBa2Cu3 06.333 ellipsoid, measured as a function of 
temperature for different values of doping starting from 
the most ordered sample with T c = 17 K. The lower 
figure describes the behavior of p s {T) very close to the 
critical temperature for two levels of doping that lead 
to the highest T c 's ( sw 17K and as 15K ) as shown in 
the upper figure. In experiments described in Ref. [l4l . 
the superfiuid density was determined from the pene- 
tration depth measurements done on the high quality 
crystals. Highly ordered samples were prepared by ex- 
tended annealing under pure oxygen gas flow, as well 
as high hydrostatic pressure at room temperature to en- 
hance CuO-chain orderin g 25 ' 26 . This procedure removes 
randomness in the spatial distribution of isolated oxygen 
vacancies that could act as pinning and scattering cen- 
ters. However, planar- type inhomogeneities are present 
as twin boundaries as these materials are not detwinned. 
In conventional superconductors, the twin boundaries are 
well known to increase locally the critical temperature^. 
In YBa2Cu3 06.333, the twin boundaries parallel to the 
c-axis are envisioned as the planes separating two re- 
gions in which the CuO chains are oriented perpendic- 
ular to each other, and are likely to serve as a source of 
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FIG. 3: (Color online). The figures show the dependence 



of bpl(T) as a function of a(T) = a(T - T c w )/T c (0) (a = 1) 
calculated from Eq. (|52p for the rescaled by means of Eq. (|55p 
parameters rii — 0.1, u m — 3W for figure a) and m = 0.1, 
u m = 5W 7 for figure b) respectively. The widths of the disor- 
der distribution W = 0.0, W = 0.1, W = 0.2 in both figures, 
are represented by the solid, dashed and dashed-dotted lines 
respectively. 



T c -enhancement as well. The reason for this is that it is 
advantageous for the oxygen vacancies to be located near 
the twin boundaries^. And the presence of narrow re- 
gions with surplus of these vacancies means the enhance- 
ment of superconductivity, and as a consequence higher 
critical temperature in the vicinity. For each plane, this 
is effectively modeled by adding the potential Eq. (|31[) to 
the conventional GL action. Randomness of u is likely 
to come from variation of the in-plane concentration of 
vacancies, stemming in its turn from the non-uniformity 
of the initial density of the twin boundaries. Since, in our 
case all such planar defects appear to increase the criti- 
cal temperature, the values of u are determined relative 
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15 20 25 30 

Temperature (K) 

FIG. 4: Measured superfluid density in the underdoped 
regime of YBa2Cu3 06.333 taken from Ref. [iH . The data 
suggest the presence of decreasing superconducting response 
above some temperature. The meaning of the curves is ex- 
plained in the text. 



to some average value, hence having both signs. Tc de- 
fined previously as the critical temperature of completely 
disorder-free sample must thus be regarded shifted up so 
that Eq. (|28p is satisfied. We should mention also that 
the issue of local superconductivity enhancement due to 
the twin-boundaries themselves has been considered pre- 
viously in Ref. [28l . The twinning planes were assumed to 
form a periodic array and be described in the GL func- 
tional by the sum of terms of the form Eq. (Q]) all having 
the same u. In our approach, we specifically highlight the 
importance of randomness in potentials for the purpose 
of the qualitative interpretation of the measured super- 
fluid density. 



It is seen from the data in Fig.|4]that a small but finite 
superfluidity persists above some temperature obtained 
by extrapolating the straight lines, describing the behav- 
ior at lower T down to zero. If we associate this value 
with Tc°^ in our approach, we can claim that the sim- 
ple model of planar disorder reproduces qualitatively the 
experimentally observed temperature dependence of p s . 
From the data, though, it is difficult to infer the precise 
value of the actual critical temperature T c . Since the aim 
of this calculation is mainly to illustrate the generic fea- 
tures resulting from randomness in T c 's, we have not at- 
tempted to determine the choices of parameters W, u m , a 
required for a precise match between the theoretical and 
experimental curves. Already this qualitative agreement 
suggests that even in experiments done on high quality 
samples, there are local regions with critical tempera- 
tures significantly above the average T c . The smallness 
of superfluid density in the tails suggests that such re- 
gions are rare and have broadly distributed local critical 
temperatures. 



IV. DISCUSSION AND CONCLUSIONS 

In this Section, we discuss the relevance of the simple 
model of random planar defects to the experimentally 
observed behavior of the superfluid density as a function 
of temperature close to T c . 

First, it is appropriate to ask what changes to the be- 
havior of p s are expected if we go to the next orders 
of expansion in defect concentration. To examine this, 
consider how one proceeds to obtain the correction to 
the superfluid density that is of second order in m . Ac- 
cording to the general strategy, it is necessary to solve 
the equation for the order parameter in the presence of 
two defects located at points Xi and X 2 - The distance 
L = \X 2 — X± \ emerges as a new parameter for the prob- 
lem, and together with the strengths of the potentials 

will determine the point at which the real positive solu- 
te 2) 

tion 'Fq '(L;x) first obtains. To find it, one can easily 
solve the corresponding generalization of Eq. for the 
Green's function Q^ 2 \X\ = 0,X 2 = L;r,r') of the oper- 
ator Lq[u] with the potential 

u(r) = —ui5(x) — U25(x — L). (55) 

The subscript (2) means now that we do all calculations 
in the presence of exactly two defects. The solution has 
the form 

G {2) (vy Z ;p,p r ) = 

5{p + P ')+u 1 C 1 {p')^u 2 C 2 {p')e-^ L 

a + 7 lyz +p 2 /(2m) ' 1 ' 



g^ { rj yz ;-p",p')^-, (57) 
-00 27T 
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C 2 (p') = J e~ ip L G {2) {Vyz;~p",p')^r. (58) 

We will not write out in full the cumbersome expression 
for the Green's function and state only that the poles of 
6 iVyz:P,p'), indicating the occurrence of a transition 
at some T > T c (0) , are determined from the equation 

. ^avassizi = . (59) 

2a 

The left hand side of Eq. (j59")l is just the corresponding 
determinant (taken at rj yz = 0) which arises in the pro- 
cess of solving the system of two coupled linear equations 
for Ci(p') and C2{p')- If wi = u 2 , Eq. (fB"9"]) reduces (upon 
the proper rescaling) to the known result^!. It is clear 
from Eq. (f5"9"|) , that for positive u\ and u 2 , the critical 
value of a = (m/2) max(iti, u^) if l-^l = °°: implying a 
complete independence of the defects. But for \L\ = 0, 
the pole occurs at a = (m/2)(ui + u 2 ) 2 , in agreement 
with the fact that if both defects are located at one point, 
their strengths simply add. It follows then that placing 
the second defect with positive u at any finite distance L, 
in addition to the one already present, increases the crit- 
ical a. This means that if there are two defects located 
not too far from each other in the sample, with strengths 
close to u m , the actual critical temperature, as a result of 
going to the second order in rii, will be higher than that 
given by Eq. (|50[) . However, this circumstance does not 
affect the qualitative interpretation of the data since, as 
was mentioned before, the presumed exponential small- 
ness of p s {T) at temperatures considerably higher than 

T c , renders the actual T c difficult to determine from 
Fig. [H Thus we will not discuss further all the calcula- 
tions to second and higher orders in rii, but rather note 
the following. Once the potential strengths are broadly 
distributed, the value of 6pj(T) at T = is not sensi- 
tive to the actual T c . To calculate it one can safely set 
u m = oo. However, the calculated b~pl(T) using Eq. ([52]) 
has a peak at a = as seen in Fig. [3J We believe that 
this non-monotonic behavior in the vicinity of Tc°\ is 



an artifact of the mean-field approximation used in our 
treatment from the very beginning. Thermal fluctuations 
when a is small, are expected to renormalize down the 
values of the superfluid density, but consideration of this 
question is beyond the scope of this paper. 

The planar defects considered in this paper, may not 
be the only ones present despite the high quality of the 
samples. Isolated and rare point defects due to oxygen 
disorder, other lattice defects such as dislocations 2 ^ are, 
in principle, not excluded^ . It is important, however, 
that among all possible types of defects, the plane-like 
ones ensure the broadest possible distribution of the local 
critical temperatures. This follows from the generalized 
Harris criterioni^i, arguing that whenever 2 — d*v > 0, 
with v being the correlation length critical exponent and 
d* the number of dimensions in which the system is ran- 
dom, disorder is relevant near the critical point. For point 
and columnar disorder, d* is equal to 3 and 2 respectively, 
whereas in the case of stacked planar defects d* = 1 . Ob- 
viously, the greater 2 — d*v is, the wider the distribution 
of local T c 's is due to randomness in parameters char- 
acterizing the defect potential^. One should remember, 
though, that not all defects, but only the extended ones 
of sizes greater than £oj can be satisfactorily accounted 
for within the framework of GL theory. Others should be 
treated using a suitable microscopic model. If necessary, 
possibly in lower quality samples, defects with spherical 
and cylindrical shapes must be considered as well, but 
these are likely to affect the superfluid density in a much 
narrower region around . Investigation of such de- 
fects is left for future work. 
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APPENDIX A: EXPANSION IN POWERS OF m 



In this Appendix, we obtain the average of the function i^jv({ui}Ri, . . . , {wjv}Rjv; r) given by Eq. (|10[) in the form 
of a series in powers of defect concentration rii. The corresponding derivation is straightforward and follows the lines 
of Ref. [l6|. We should notice first, however, that the function in angular brackets in Eq. pU|) is not symmetric, but 
the final result will not change if we replace Fn({ui}Ri, . . . , {uyv}Riv; r) with the function 

■^({-"l}, • • ■ , {u N }; R-i, R»; r) = -^P [{wi}, . . . , {u N }} Fjv({ui}Ri, . . . , {uatIRtv; r), (Al) 



13 



symmetrized over all sets {ui} for a particular location of the defect. ( P is the symmetrization operator) This creates 
functions symmetric with respect to permutations of coordinates Hi , . . . , Rjv ■ It is possible to verify then that for 
any finite number of such functions F sm ({ui}, . . . , {u m }; Ri, . . . , R m ; r), (m > 0) 

F sA r({ui}, . . . , {u N }; Ri, , Rat; r) = F (r) + ^ $i{{u t }; R,; r) 

i 

+ J2 [$ s2 ({M> {%•}; R, Rj ■; r) - $i(KI; r) - $i(W; R,; r)] 

i<j 

+ J! [ $ s3({wi},{%-},{u fc };Ri,Rj,R fe ;r) - $ s2 ({ui},{%};Ri,Rj;r) - $ s2 ({wJ, {w fc }; R l; R fc ; r) 

- $ s2 ({%}, {ufc}; Rj, R/t; r) + R 4 ; r) + $i({^}; Rj; r) + $i{{u k }; R fc ; r)] + • • • (A2) 

In Eq. flMJ (to > 2), 

$ S m(W, ■ ■ ■ , {u m };Ri, . . . ,R m ;r) = F sm ({ui}, {u m };Ri, . . . ,R m ;r) - F Q (r), (A3) 

and these $ S m({Mi}, . . . , {u m }; Ri, . . . , R m ; r) are also symmetric. Fq(t) is a function of r only, and it, as well as 
$i({ui}; R^ r), does not require symmetrization. Assuming that all defects are located at points Ri, and that Fq(t) 
is the value in the absence of defects, we can take the thermodynamic limit N — > oo, V^-^oo, n, = N/V in Eq. (|A2j) 
to get 

(^({uijRi, • • • , {u n }Rn; r)> = F (r) + n, J $i({iti}; Ri; r)dRi 
+ // [ $s2 ^ Ml ^' ^ U2 ^ Rl ' R2; r ) - R ii r ) - *i({ M 2}; R 2 ; r)] w 2 ({^i}, {u 2 }; R 2 - Ri) dRidR 2 

+|y- [$ S 3({wi}, {" 2 }, {w 3 }; Ri, R 2 , R 3 ; r) - $ s2 ({ui}, {ua}; Ri, R 2 ; r) - $ s2 ({w 2 }, {u 3 }; R 2 , R 3 ; r) 

-$ s2 ({u 3 }, {^i}; R 3 , Ri; r) + Ri; r) + $i({u 2 }; R 2 ; r) + $i({w 3 }; R 3 ; r)] • 

w 3 ({wi}, {w 2 }, {m 3 }; R 2 - Ri, R 3 - R 2 ) dRidR 2 dR 3 + • • • (A4) 

The full average then follows in a straightforward way by substituting this average over positions into Eq. (jlOp . In 
Eq. (|A4|) . u> 2 ({ui}, {u 2 }; R 2 — Ri) is the probability of the simultaneous occurrence of two defects at points Ri and R 2 
characterized by the sets {ui} and {u 2 } respectively. Obviously, iu 2 ({ui}, {it 2 }; R 2 — Ri) = u> 2 ({iti}, {m 2 }; Ri — R 2 ). 
Probabilities 11)3, 11)4, . . . in higher order terms have analogous meaning, and in general depend on {ui}. However, 
this dependence is likely to be noticeable only if defects are close to each other. All functions w m will depend on 
concentration rii, but in the limit all |Rj — R/-| 3> l/(rii) 1 ^ d , they quickly tend to unity. Those functions must also 
factorize if one of the coordinates tends to infinity, meaning, for example, that ui 3 ({ui}, {u 2 }, {w 3 }; R 2 — Ri,R 3 — 
R 2 ) -> uj 2 ({ui},{u 3 };R 3 - Ri), if |R 2 | -> 00. 
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